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ABSTRACT 

We measure the number of companions per g alaxy (N r ) as a funct ion of r-band absolute magnitude 
for both the Sloan Digital Sky Survey and the ICroton et al.l (|2006l ) semi-analytic catalog applied to 
the Millennium Run simulation. For close pairs with projected separations of 5-20 h~ x kpc, velocity 
differences less than 500 km s _1 , and luminosity ratios between 1:2 and 2:1, we find good agreement 
between the observations and simulations, with iV c consistently close to 0.02 over the range —22 < 
M r < —18. For larger pair separations, N c (M r ) instead becomes increasingly steep towards the faint 
end, implying that luminosity-dependent clustering plays an important role on small scales. Using the 
simulations to assess and correct for projection effects, we infer that the real-space N c (M r ) for close 
pairs peaks at about M*, and declines by at least a factor of two as M r becomes fainter. Conversely, by 
measuring the number density of close companions, we estimate that at least 90% of all major mergers 
occur between galaxies which are fainter than L* . Finally, measurements of the luminosity density 
of close companions indicate that L* galaxies likely dominate in terms of the overall importance of 
major mergers in the evolution of galaxy populations at low redshift. 

Subject headings: galaxies: evolution, galaxies: interactions, surveys, galaxies: statistics 



1. INTRODUCTION 

Galaxy mergers can produce dramatic changes in the 
morphological, nuclear and star forming properties of 
galaxies over relatively short timespans. As a result, 
mergers have been invoked to explain a number of as- 
pects of galaxy evolution. In recent years, large redshift 
surveys have paved the way for systematic observational 
studies of candidate mergers, while theoretical modelling 
of structure formation has produced key advances in our 
understanding of the role of merging in a cosmological 
context. In general, these efforts have focussed on two 
key aspects of merging: (1) the effects of merging on the 
constituent galaxies, and (2) the frequency with which 
merging occurs, as described by the merger rate and re- 
lated quantities. 

An increasingly popular method of identifying can- 
didate merging systems is the use of close galaxy 
pairs, which are the precursors to mergers. With 
careful choices of close pair criteria, and correction for 
projection effects (ie., contamination by non- merging 
pairs), the frequency of close pairs should correlate 
with the merger rate. Recent studies using cosmological 
simulations support this idea, demonstrating that 
most close pairs merge on relatively short timescales 
(|Kitzbichler fe White! 120081) . In addition, the prop- 
erties of paired galaxies can provide insight into the 
nature of merging galaxies both before and during 
the encounter. Galaxies in close pairs have higher 
asymmetri es than galaxies in wider separation pairs or 
the field (iHernandez-Toledo et ail 120051 : IPatton et all 
l2005t IDe Propris et all 120071 ). confirming that interac- 
tions and mergers are prevalent in these systems. Star 
formation is enhanc e d in close pairs a t low redshift 
(iLambas et al.l 120031: lAlonso et all l2004t iNikolic et al l 
12004 IPatton et all 120051 : iSol Alonso et alj 120061 : 
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been triggered by mergers or interactions. Differences in 
metallicities between paired and field galaxies are consis- 
tent with a scenario in which interactions funnel gas to 
the centr al regions of galaxies involved in these c lose en- 
counters (jKewlev et all 120061; lEllison et alj 120081 ). Most 
ultraluminous infrared galaxies (ULIR Gs) originate 
from m ajor mergers of gas rich galaxies (|Dasvra et al.l 
2006alrbr ). while approximately half of the luminous 
infrared galaxies (LIRGs) at low r edshift appear to be 
undergoing interactions or mergers dWang et al. ll2006). 

Using close galaxy pairs and/or galaxy asymmetries 
as indicators of imminent or recent mergers, the merger 
rate and its evolution has now be en measured using a 
number of large redshift surveys JCa rlbcrg et alJl2000t 
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Evolution estimates range from roughly (1 + z) 5 to 
(1 + z) 3 , implying widely differing scenarios at z ~ 1 
and above. At least some of these discrepancies result 
from the use of different pair criteria. For example, sim- 
ulations and semi-analytical models of galaxy formation 
indicate that the merger rate and its evolution depends 
on factors such as galaxy mass, pair mass ratio and envi 



ronmcnt (Khochfar & Burkcrt 2001; Bcrricr ct al. 



2006 



2008 



Mailer et a l. 2006; Co x et all 120081: iGuo fc Whit 
Kitzbichler fc Whitd 120081 ). 

In order to better understand the role of merging, 
and to reconcile merger rate measurements from dis- 
parate samples, we need to assess how the frequency 
and nature of merging depends on factors such as en- 
vironment, mass ratio (ie., major versus minor mergers), 
properties of the progenitor galaxies (e.g., dry mergers 
versus gas-rich mergers), and the mass (or luminosity) 
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of the merging galaxies or merger remnants (e.g., for- 
mation of massive galaxies versus L* galaxies). Sig- 
nificant observational progress has been made in all 
of these areas in recent years. Galaxy gr oups appear 
to be an ideal environment for mergers rtGotol [20051: 
Brough et all 120061: iMiles et al l [200llRobotham et al 



2006J iZandivarez e"t all I2006t ICoziol fc Plauchu-Frav . 
2007t iNolan et all 120071). though the infall regions of 
rs (|va 



clust ers( |van Dokk um et alJll99a |Tran et al.112 005: Moss 



120061) and the low density field (jBarton et al J 120071) are 
important too. Induced star formation appears to be 
strongest in major mergers, or in t he lower luminosity 
(or mass) members of minor mergers dWoods et al.l f2006 : 
iWoods fc GelTerl [20071 : lEllison et alJ |2OO0). Dry merg- 
ers have been invoked to explain th e assembly of mas- 
sive elliptical ga l axies since z ~ 1 (|van Dokkuml 12005: 
iBell et al.|[2006bl:lNaab et al.ll2006D . though this process 
likel y cannot explain a l l recently formed ea r ly type galax- 
ies ^ Coxet "all 120061: fBrown et al.l 120071 : feuridy et all 
120071 : IScarlata et al.l[200l . 

In this paper, we investigate the luminosity depen- 
dence of the merger rate at low redshift, using close 
galax y pairs in the S loan Digital S k y Surv ey ijYork et all 
2000), and in the ICroton et al.l (|2006l ) semi-analytic 
galaxy catalogs derived from the Millennium Run sim- 
ulation ([Springel et al.ll2005[ ). Both of these samples are 
large enough that, in addition to being able to measure 
close pair statistics as a function of absolute magnitude, 
we also have the luxury of being very selective in how we 
choose our close pairs. In particular, we require all of our 
pairs to have spectroscopic redshifts for both members, 
projected separations less than 20 h~ x kpc, and relative 
velocities less than 500 km s _1 . In addition, we require 
our companion sample to be volume limited for all lu- 
minosity ratios between 1:2 and 2:1, thereby providing a 
cleaner match to the major merger candidates we seek 
to identify. Moreover, our measurements are carried out 
in the r-band, yielding absolute magnitudes that are a 
better proxy for stellar mass than those at shorter wave- 
lengths. This is beneficial for close pair studies, since 
merger-induced star formation is likely to affect the lumi- 
nosities of galaxies in pairs more than normal (isolated) 
galaxies. 

We begin by describing the creation of our SDSS spec- 
troscopic and photometric samples in §0 Section [3] out- 
lines the calculation of the number of close companions 
per galaxy (N c ) for SDSS, including corrections for spec- 
troscopic incompleteness. We make a direct comparison 
with the Millennium Run simulation in § [4[ and derive 
real space pair statistics for both SDSS and Millennium 
in §[5[ With additional assumptions, we then relate these 
close pair statistics to the merger rate in § [7J] We sum- 
marize our conclusions in § [71 Throughout this study, we 
adopt cosmological parameters of fl m = 0.3, f^A = 0.7, 
and Hq = 100 h km s _1 Mpc -1 . For brevity, we express 
all absolute magnitudes as M r instead of M r — 51og(/i). 

2. DATA 

The fifth data release (DR5) of the Sloan Digital Sky 
Survey (hereafter SDSS) consists of ugriz imaging span- 
ning over 8000 deg 2 , along with spectra of about one 
million galaxies, quasars, and stars within a 5713 deg 2 
subse t of the imaging area (|Adelman-McCarthv et all 
2007). In this study, we use the main galaxy sample, as 



described bv lStrauss et al.l (I2002D. Un l ike th e luminous 
red galaxy sample of lEisenstein et al.l (|2001| ). this sam- 
ple is designed to be independent of galaxy luminosity 
and Hubble type, and therefore probes a representative 
sample of galaxies, including both gas-rich and gas-poor 
( "dry" ) merger candidates. 

Our goal is to carry out a census of galaxies with close 
physical companions. We restrict our observational sam- 
ple to close galaxy pairs for which redshifts are available 
for both galaxies. This reduces contamination due to 
non-merging systems, and allows one to measure rest- 
frame properties of the pairs and their member galaxies. 
However, the SDSS is not a complete spectroscopic sam- 
ple, and in fact the minimum fiber separation of 55" bi- 
ases the sample against the close pairs we are interested 
in. Fortunately, it is possible to measure and correct for 
spectroscopic incompleteness by using the photometry 
of all galaxies, regardless of whether or not spectra were 
obtained. To this end, we now describe the creation of 
both spectroscopic and photometric catalogs of galaxies 
drawn from the main galaxy sample. 

2.1. Spectroscopic Catalog 

The SDSS main galaxy sample consists of galax- 
ies with r-band Pe trosian magnitudes of m r < 17.77 
(jStrauss et alj [2002) , after correction for Galactic extinc- 
tion. However, regions covered by the first data release 
of the survey had limi ting apparent magnitu des rang- 
ing from 17.5 to 17.77 ([Abazaiian et alj |2003). In addi- 
tion, the sample becomes incomplete at the bright end, 
where aut omated deblending of large galaxies becomes 
unreliable ([Strauss et al.ll2002l ). and introduces concerns 
about single galaxies mistakenly being classified as close 
pairs. To ensure a complete and reliable sample, we 
therefore begin by restricting our catalogs to the range 
14.5 < m r < 17.5. We also ensure that every galaxy has 
a measured spectroscopic redshift; moreover, we require 
the SDSS "zConf" parameter to be at least 0.7, thereby 
ensuring that the confidence in each redshift measure- 
ment is at least 70% (in most cases, it is much higher). 
Using these criteria, we create a preliminary spectro- 
scopic sample by querying the SDSS online database 2 . 

For every galaxy in the spectroscopic sample, we mea- 
sure the r-band absolute magnitude at redshift z, as 
given by 

M r = m r - 51ogd L (z) - 25 - k r - E(z), (1) 

where m r is the extinction-corrected Petrosian r-band 
magnitude, oIl is the luminosity distance, and k r and 
E(z) are the fc-corrections and passive stellar evolution 
corrections, respectively. We measure fc-corrections using 
the SDSS ugriz photome try, employing the kcorrec t soft- 
ware ( version vl-1-4 ) of IBlanton fc Roweid (|2007f ). Fol- 
lowing [PattoK^iE^n (|2002f ) . we parameterize the evolu- 
tion correction as E(z) = —Qz, where Q is determined 
from measurements of the galaxy luminos ity function 
(LF) and z is the redshift (|Lin et al.|[l999l ). We adopt 
Q = 1.8 as a n average of the SPSS 01 r and 01 g Q mea- 
surements of IBlanton et alJ (|2003al ) , since rest-frame r 
lies between these two passban ds (this agrees well wi th 
the evolutionary correction of iTegmark et alJ (|2004h ). 
Given the relatively low redshift of our sample, however, 

2 http:/ /casjobs. sdss.org/CasJobs/ 
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this evolution correction is small, and removing it en- 
tirely (ie., Q = 0) does not significantly change any of 
the conclusions in this paper. 

2.2. Photometric Catalog 

In order to measure and correct for incompleteness in 
our spectroscopic catalog, we also create a photometric 
catalog, in which galaxies satisfy the same flux limits as 
the spectroscopic catalog, but are not required to have a 
measured redshift. The sky area covered by DR5 imag- 
ing is abo ut 40% larger than the area with spectroscopic 
coverage (jAdelman- McCarthy et al.ll2007[ k to match the 
footprints of our spectroscopic and photometric catalogs, 
we therefore also require every galaxy in both samples to 
have at least one galaxy with a spectrum within 12 ar- 
cminutes of its position (excluding itself). 

2.3. Spectroscopic Incompleteness 

By comparing our resulting photometric and spec- 
troscopic catalogs, we find that the average spectro- 
scopic completeness is 88%. While spectroscopic tar- 
get selection was desi gned to provide uniform coverage 
(jBlanton et al.ll2003bl ). the completeness does vary con- 
siderably from one part of the sky to another, with the 
completeness in some regions falling well below the mean, 
while in other regions (particularly those covered by more 
than on e SPSS plate ), the c ompleteness is considerably 
higher. iPatton et alJ (|2002h demonstrated that the ob- 
served number of companions per galaxy scales with the 
spectroscopic completeness. While this bias can be cor- 
rected for, one can also minimize its impact by excluding 
regions with low spectroscopic completeness. Such an ex- 
clusion will also remove galaxies which lie close enough 
to the survey boundaries that close companions will be 
missed. With this in mind, and with a desire to mea- 
sure and correct for the remaining spectroscopic incom- 
pleteness, we compute a measure of local spectroscopic 
completeness for every galaxy in our spectroscopic and 
photometric catalogs. We consider an area around each 
galaxy within an outer radius of one degree, and an in- 
ner radius of 55" (the latter corresponds to the minimum 
fiber separation) . After counting the number of enclosed 
galaxies in the spectroscopic and photometric samples, 
we take the ratio of these two numbers, denoting this 
quantity f s . We then require every galaxy in our cata- 
logs to have f s > 0.75. This restriction excludes 4.7% 
of the galaxies in the spectroscopic catalog. The main 
conclusions of this paper are unchanged if we instead use 
f s > 0.7 or f a > 0.8 (excluding 2.5% and 9.5% of the 
spectroscopic catalog respectively). 

3. SDSS CLOSE PAIR STATISTICS 

3.1. Methodology 

We now proceed to identify close galaxy pairs from 
our spectroscopic sample. We measure three key quan- 
tities for each galaxy pair: projected physical separation 
r p , rest- frame relative velocity At>, and absolut e mag- 
nitud e difference AM r . Following fPatton et ail (|2000l 
2002), we define a "close pair" to have 5 < r p < 20 h 
kpc and Av < 500 km s -1 . On order half of the pairs 
satisfying these criteria are known to exhibit morpho- 
logic al signs of interacti ons, based on visual classifica- 
tion (|Patton et al.l l2000f) and quantitative measures of 



asymmetry (jPatton et all 120051 : |Pe Propris et al.ll2007f) . 

In addition, in order to focus on major merger candi- 
dates, we require that |AM r < 0.753, ensuring a pair 
luminosity ratio between 1:2 and 2:1. This criterion is 
preferable to the more common approach of selecting 
both host and companion galaxies from the sam e fixed - 
range in absolut e mag n itude (e.g., IPatton et al.1 (|200Ct 
120021) : iLin et alJ (|2004f ): iDe Propris et alJ (|2005[ )). since 
that approach includes a wider range in luminosity ra- 
tios, and becomes increasingly incomplete at luminosity 
ratios significantly different from 1:1. 

3.2. Choosing Potential Host and Companion Galaxies 

Figure [T] contains a plot of absolute magnitude ver- 
sus redshift for 10,000 galaxies selected at random from 
our spectroscopic sample of about 337,000 galaxies. We 
begin our sample selection by measuring the bright and 
faint limits in absolute magnitude as a function of red- 
shift within which galaxies of all spectral types will have 
14.5 < m r < 17.5. These are shown by the upper and 
lower dashed (black) curves in Figure [TJ and are com- 
puted using estimates of the minimum and maximum 
k-corrections respectively. We denote all galaxies lying 
within these limits as potential companion galaxies. 

Within this sample, we then identify the subset of 
galaxies for which all companions (with \ AM r \ < 0.753) 
are detectable. These galaxies, which we refer to as po- 
tential host galaxies, are contained within the two solid 
(red) curves in Figure[T] Aside from spectroscopic incom- 
pleteness (which we correct for statistically), this pro- 
vides a sample of host galaxies for which we can carry out 
a volume limited search for companions. It follows that 
our sample will contain close pairs in which either one or 
both galaxies fall into the sample of potential hosts. 

We wish to compute pair statistics as a function of M r , 
over as large a range as feasible. However, in order to en- 
sure completeness, we must impose some restrictions on 
the range of host galaxy luminosities. At the bright end, 
Figure[T]indicates that our sample contains relatively few 
galaxies close to M r ~ —23. Given the need to detect 
companions which are 0.753 magnitudes brighter than 
any potential host galaxy, we therefore impose a mini- 
mum M r = —22. At the faint end, the key issue is to 
decide on a reasonable minimum redshift to use, as this 
will dictate a maximum allowable absolute magnitude. 
Galaxies at the lowest redshifts have the least certain 
absolute magnitudes, due to the increased influence of 
peculiar velocities on the observed redshift. In addition, 
such nearby galaxies are typically the most challenging 
to obtain accurate photometry for, since the SDSS de- 
blending algorithm tends to break down more often for 
galaxies with large apparent sizes. With these consider- 
ations in mind, we impose a conservative maximum M r 
of —18 on our sample of potential host galaxies, corre- 
sponding to a minimum redshift of 0.015. Companions 
are permitted to lie at slightly lower redshifts due to the 
allowed velocity difference of 500 km s" 1 . 

With these criteria, we find a total of 477 host galax- 
ies which have at least one close companion. In order 
to ensure that the SDSS pipeline has been successful in 
detecting real galaxy pairs, we visually inspect the SDSS 
images of every detected host galaxy. The contamination 
was found to be very small, with only 0.8% of the host 
galaxies being spurious. The affected systems consist of 
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one edge-on disk, and one ring galaxy: in both cases, 
the SDSS pipeline mistakenly detected two galaxies. We 
therefore remove these galaxies from our sample, leaving 
473 host galaxies in the range —22 < M r < —18. The 
basic properties of these galaxies are listed in Table [1] 
(this table contains a subset of the table, which is to be 
published in its entirety in the electronic version of the 
journal). 

A significant fraction of these galaxies exhibit morpho- 
logical signs of interactions, though we defer a more rig- 
orous structural analysis to a future paper. We note 
that this sample is more than twice as large as the Mil- 
lenniu m Galaxy Catalogue sample of iDe Propris et al.l 
( 2005) , despite our more rigid requirement that the sam- 
ple be volume-limited for close companions which are 
major merger candidates (\AM r < 0.753|). 

3.3. Small Scale Spectroscopic Incompleteness 

In Section [2. 3[ we described our algorithm for measur- 
ing the local spectroscopic completeness for every galaxy. 
This provides a measure of completeness on scales of on 
order half a degree. However, constraints on fiber place- 
ment require angular separations of at leas t 55" between 
any t wo targets assigned to the same plate (|Strauss et al.l 
2002). As a result, spectroscopic completeness drops 
sharply at smaller angular pair separations. At a redshift 
of 0.1, 55" corresponds to a projected separation of 71 
/i _1 kpc, meaning that most of the close pairs of interest 
in this study (r p < 20 hr 1 kpc) are found at these small 
angular separations. Fortunately, most plates contain re- 
gions of overlap with adjacent plates, and some regions 
are observed using two or more plates. As a result, the 
spectroscopic sample contains enough close angular pairs 
that we are able to model the i ncompleteness and c orrect 
for it, following the method of iPatton et al.l (|2002f ). 

First, we measure the ratio of spectroscopic to photo- 
metric pairs (N zz /N pp ) as a function of angular separa- 
tion 9, as shown in the upper panel of Figure O For a 
fair sample, in which spectroscopic completeness is in- 
dependent of pair separation, one would expect to find 
N zz /Npp ~ fg, where f s is the overall spectroscopic com- 
pleteness of the survey. Instead, we see a sharp drop 
in N zz /Np P below 55", as expected. We model this in- 
completeness by fitting a function g{9) to these data. 
We then multiply each spectroscopic pair by a weight 
wg = fg/ g{9) . The resulting corrected values of N zz /N pp 
are plotted in the lower panel of Figure [2 It is clear 
that this weighting scheme is successful in removing the 
angular dependence of the small scale spectroscopic in- 
completeness, with the corrected N zz /N pp ks /| at all 
relevant angular separations, to within the reported er- 
rors. Therefore, by applying wg weights to each detected 
pair, we can remove this very significant bias from our 
measurements . 

3.4. N c (M r ) for SDSS Pairs 

We now have a catalog of host galaxies for which we 
can measure close pair statistics, using statistical weights 
to correct for spectroscopic incompleteness . Weights are 
combi ned using the approach outlined by IPatton et al.l 
(2002). First, the number of companions for galaxy i, 
summed over all companions j, is given by 

•v. Y^L '«•„ • (2) 

j 



The statistical weights in this equation are used to cor- 
rect the observed number of companions to the number 
that would have been observed in a complete redshift 
survey. The mean number of companions per galaxy, 
weighted by spectroscopic completeness, is then given by 

N c = E i/ 5 - 7Vc ' . (3) 

This weighting scheme places greater importance on 
galaxies in regions of higher spectroscopic completeness, 
thereby minimizing statistical uncertainties. 

The resulting measurements of N c (M r ), calculated us- 
ing equations and [3J are indicated by the solid (blue) 
curve in Figure [3] In addition, to facilitate future com- 
parisons with other surveys, we provide these results in 
the first row of Table [2] We find that N c (M r ) is ap- 
proximately constant over the given range in absolute 
magnitude, with a mean of N c = 0.021 ± 0.001. Given 
that few galaxies have more than one close companion, 
this implies that ~ 2% of galaxies with —22 < M r < — 18 
have a close companion, independent of M r . 

Overall, our mean value of N c agrees quite well with re- 
lated measurements in the literature. For example, using 
identi cal r v and Ay criter ia, but no limits on luminosity 
ratio. IPatton et alj (|2000fl find iV e (-21 < M B < -18) = 
0.0226 ± 0.0052 at z = 0.015, while IDe Propris et all 
(p0"05h find N c = 0.0357±0.0027at z = 0.116. Given our 
additional requirement that companions be within a fac- 
tor of two in luminosity of their host galaxy, both results 
are broadly consistent with the somewhat lower N c that 
we find. Our mean N c appears to be several time s large r 
than the SDSS pair fraction of iKartaltepe et alj (2007). 
However, their close pair criteria are substantially dif- 
ferent from ours; in particular, they require My > —20 
for both members, they do not require spectroscopic red- 
shifts for both me mbers of their pairs, and their sample 
is derived from the lAllam et al.l (|2004l ) sample of merg- 
ing pairs, which have a number of additional criteria im- 
pose d (including isolation). A meaningful comparison 
with IKartaltepe et alj (|2007f ) is therefore not feasible. 

The flat shape of N c (M r ) that we find is a surpris- 
ing result. Based on the LF alone, one would expect to 
find more companions (of comparable luminosity) close 
to intrinsically faint galaxies, since dwarf galaxies are 
much more numerous than giants. In other words, N c 
shoul d be proportional to the number density of galax- 
ies (|Berrier et aJ J 12006) . However, the number of com- 
panions per galaxy is s ensitive to both number density 
and c lustering strength ([Patton et al.ll2000l : lBerrier et al.l 
2006). Cl ustering strength is known to increase with lu- 
minosity (iNorberg et al1l200lL 120021 : IZehavi et aDl2002t 
iLi et al.ll2006h . This provides a competing effect, which 
must be comparable in size to the density effect if it is to 
explain the flat N c (M r ) that we find. However, at sep- 
arations of 1 00 h~ l kpc, th e correlation function mea- 
surements of ILi et all (|2006f ) indicate that the clustering 
strength of SDSS galaxies is roughly independent of ab- 
solute magnitude for M r > —21, and is approximately 
twice as high in the range —22 < M r < —21. In order to 
explain the roughly flat N c (M r ) we detect over the range 
—22 < M r < —18, the luminosity dependence of cluster- 
ing must be considerably stronger on the smaller scales 
relevant for our close pairs (~ 10 hr 1 kpc). However, we 
must also rule out any potential luminosity-dependent 
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biases in our sample, which could in principle contribute. 

It is possible to break the degeneracy between den- 
sity and clustering effects by considering wider separa- 
tion pairs, since clustering strength diminishes with pair 
separation. We therefore compute N c (M r ) as a func- 
tion of pair separation out to 100 hr 1 kpc (see Figure |3] 
and Table [J). As expected, N c (M r ) steepens towards 
the faint end as pair separation increases. This implies 
that N c (M r ) scales with number density on large scales, 
but flattens out on small scales as a result of the added 
effects of luminosity-dependent clustering. We conclude 
that the roughly flat N c (M r ) that we observe for close 
pairs is seen only on the small scales relevant for galaxy 
interactions and mergers, and is unlikely to be due to any 
overall luminosity-dependent bias within our sample. 

4. MILLENNIUM CLOSE PAIR STATISTICS 

Given the somewhat surprising trends seen in our 
SDSS pair statistics, and to facilitate comparison with 
theoretical models of galaxy formation and evolution, 
we now apply our techniques to a sample drawn from 
the semi-analytic galaxy catalogs of lCroton et al.l {2006) , 
which were created using the output of the Millen- 
nium Run simulation (|Springel et al.ll2005f ). These cat- 
alogs have been shown to reproduce the overall proper- 
ties of galaxies in the local universe, including the lu- 
minosity function, the two point c orrelation function , 
and the pairw ise velocity dispersion (jCroton et al. I [20061 
iLi et al.ll2007l ). The simulations provide the additional 
benefit of three dimensional positions and velocities, 
which we will use to assess and remove projection ef- 
fects. While the reader is referred to ICroton et al.l 
(2006) and references therein for a detailed description 
of these catalogs, we note that merging is treated by fol- 
lowing dark matter subhalos down to a mass limit of 
1.7 x 10 10 h~ 1 M(^, and then using the d ynamical friction 
formula of lBinnev and Tremaind l) 19871 ) to estimate the 
remaining time until a merge r takes place. 

For this analysis, we use the lCroton et al.1 (|2006f ) ugriz 
catalog, which is complete for galaxies more luminous 
than M r — —16.6. The stated resolution is 5 /i -1 kpc, 
which conveniently coincides with our minimum pair sep- 
aration criterion. We confirm, however, that the two- 
point galaxy correlation function has the expected power 
law form to below this level, so there should not be any 
unexpected effects near the resolution limit. 

4.1. Mock Redshift Catalogs 

The ICroton et al.l (|2006| ) catalogs contains three di- 
mensional positions and velocities for approximately 9 
million galaxies at redshift zero, within a cube which is 
500 hr 1 Mpc on a side. In order to make a direct com- 
parison between Millennium and SDSS, it is necessary 
to transform this real-space catalog into a mock redshift 
catalog. To this end, we begin by placing the observer at 
the origin of the Millennium cube, and then computing 
the corresponding right ascension, declination, redshift, 
and apparent magnitude of each galaxy in the sample. 
In order to match our observed SDSS sample, we then 
apply the same flux limits (14.5 < m r < 17.5). In addi- 
tion, we impose a minimum redshift of 0.022, to ensure 
that we sample only galaxies which are more luminous 
than M r = —16.6 (the completeness limit of the Cro- 
ton catalog). Finally, we impose a maximum redshift 



of 0.17, which ensures that we do not probe distances 
greater than 500 hr 1 Mpc (the size of the simulation 
cube). The resulting catalog contains ro ughly 300,000 
galaxi es, or about one thirtieth of the full ICroton et al.l 
(2006) sample. 

There are two disadvantages to this simple approach. 
First, the vast majority of the galaxies in the simulation 
are discarded as a result of the imposed flux limits, de- 
spite the fact that they are sufficiently luminous to be 
of interest. Secondly, in order to assess projection ef- 
fects, it would be useful to view galaxy associations from 
a variety of vantage points. 

Fortunately, the periodic boundary conditions imposed 
on the simulations provide a way forward. Follow- 
ing part of t h e "ra ndom tiling" technique outlined by 
iBlaizot et al.l (|2005t ). we generate a suite of mock red- 
shift catalogs. Specifically, for each mock catalog, we be- 
gin by applying to the Millennium cube a translation in 
each of the three spatial directions, with the size of each 
translation being a random fraction of the box size (500 
h~ 1 Mpc). For any resulting coordinate value greater 
than the box size, we subtract the box size, thereby 
"wrapping around" the cube. Finally, we rotate the cube 
around each of the three spatial axes by a random (in- 
teger) multiple of 7r/2. Using this approach, we create 
a suite of 30 mock redshift catalogs, sampling a total of 
about 9 million galaxies (comparable to the number of 
galaxies in the original data cube). 

4.2. Redshift- space Pair Statistics 

We then proceed to compute pair statistics on all of 
the mock redshift catalogs, selecting host and compan- 
ion galaxies in the same manner outlined in Sections 13. II 
and 13.21 One key difference in technique must be noted, 
however: given that the Millennium data is a redshift 
zero realization, we compute co-moving rather than phys- 
ical projected separations in order to recover the correct 
z = physical separations. 

Our results are presented in Figure [4] and in the first 
two columns of Table |U Overall, the trends seen are 
broadly similar to those for SDSS shown in Figure [3) 
The distribution is again relatively flat for the small- 
est separation pairs, and increases towards the faint end 
as pair separation increases. This implies that the pro- 
jected two point correlation becomes steeper on small 
scales as luminosity increases, which is consistent with 
the increase in steepn ess with stellar mass f o und f or Mil- 
lennium galaxies by iKitzbichler fc White! (|2008f) . For 
close pairs (5 < r p < 20 h^ 1 kpc), N c {M r ) peaks at 
a value of about 0.02 at M r ~ —20.75 . Given that 
M* is approximately —20.6 in the r filter (Blanton et al. 
l2003al) . this implies that M* galaxies are the most likely 
to have close companions. For close pairs over the range 
-22 < M r < -18, the mean N c is 0.0183±0.0001, which 
is about 15% smaller than we found for SDSS. The most 
significant difference between SDSS and Millennium is 
seen for the larger separation pairs, in that the Millen- 
nium N c (M r ) is considerably steeper. We discuss a pos- 
sible explanation for this in Section 16.11 However, given 
that we are primarily interested in the smallest sepa- 
ration pairs, the general agreement seen between SDSS 
and Millennium in this regime justifies further compari- 
son between these samples. 
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4.3. Projection Effects 

One of the most challenging aspects of close pair stud- 
ies is the contamination of pair samples due to projection 
effects. For any observed galaxy pair, even when spec- 
troscopic redshifts are available for both galaxies, one 
cannot be certain that the galaxies are close enough to 
merge. In addition, without measurements of transverse 
velocities, it is also possible that the relative velocities of 
the member galaxies are too high for coalescence to oc- 
cur. For pairs with 5 < r p < 20 h~ x kpc and Av < 500 
km s -1 , iPatton et alj (|2000f) estimated this contamina- 
tion to be on the order of 50%. One must apply such 
a statistical correction for projection effects in order to 
relate measured pair statistics to their real-space (three- 
dimensional) equivalents. 

The three dimen sional information available in the 
ICroton et al.l (|2006f ) catalogs enables us to instead mea- 
sure the level of contamination directly. For every pair 
of galaxies that we observe in our mock redshift catalogs, 
we can measure the three dimensional physical separa- 
tion, along with the three dimensional relative velocity. 
If these quantities are less than 20 h^ 1 kpc and 500 km 
s -1 respectively, it is likely that a merger is imminent 
(jPatton et alj 12000). In addition, every pair satisfying 
these three dimensional criteria will also satisfy the red- 
shift space criteria (i.e., in terms of r p and Av) when 
viewed from any vantage point. 

For our Millennium mock catalogs, we measure the 
fraction of close companion s which satisfy the t hree di- 
mensional criteria; following IPatton et al.l (|2000D . we re- 
fer to this fraction as f^o- Over the range —22 < M r < 
-18, we find an overall f 3D of 47.4% ± 0.3%. This com- 
pares well w i th th e rough estimate of fao ~ 50% of 
IPatton et al.l (2000), although their estimate applies to 
different absolute magnitude criteria (— 21 < Mb < — 18 
for both hosts and companions). 

In addition to measuring the overall /sd, it is also use- 
ful to probe the dependence of this quantity on absolute 
magnitude. If it varies significantly with M r , this will 
have implications for any merger rate estimates we wish 
to glean from our N c (M r ) measurements. The results are 
given in Figure El and in the third column of Table |U In- 
terestingly, /3£> is seen to vary by a factor of about three 
over the range in M r probed, decreasing from ~ 56% 
at M r < M* to about 17% at M r 18. This im- 
plies that most of the lowest luminosity close pairs in 
our Millennium mock catalogs are not candidates for im- 
minent mergers, while at least half of the close pairs with 
M r < M* are likely to undergo imminent mergers. 

5. CLOSE PAIR STATISTICS IN REAL SPACE 

5.1. Millennium 

In the preceding section, we reported our measurement 
of fzD{M r ) from the Millennium simulation. Multiply- 
ing this function by the redshift space close pair statistics 
reported in ^ 14.21 we arrive at real space pair statistics, in 
which projection effects have been removed. In Figure [6) 
we present these results for the Millennium simulation, 
along with the redshift-space results given earlier. We 
also provide tabulated values in column 4 of Table |U In 
real space, we find that N c {M r ) peaks at M r ~ M*, at 
a value of ~ 0.011, rather than ~ 0.02 in redshift space. 
Also, rather than being relatively flat (as it is in redshift 



space), we now see a strong decline towards fainter lu- 
minosities, as a result of the associated decline in f 3 n- 
The real space N c drops to about 0.002 at the faintest lu- 
minosities probed, which is much lower than in redshift 
space. These results imply that L* galaxies are much 
more likely to have close companions than 0.1L* galax- 
ies. This striking difference is not apparent in redshift 
space, due to the presence of projection effects. We also 
see a hint that the real space N c declines at the bright 
end (M r ~ —23), though this is seen only in the most 
luminous bin in the simulations. 

5.2. SDSS 

The strong contribution from projection effects seen 
in the Millennium simulation has important implications 
for the SDSS redshift space pair statistics presented in 
§ 13.41 In particular, the relatively flat N c (M r ) observed 
for SDSS close pairs almost certainly does not hold true 
in real space. Without three dimensional positions and 
velocities for SDSS galaxies, we are unable to directly 
remove projection effects. However, given the general 
agreement between the SDSS and Millennium redshift 
space pair statistics (see § 14. 2|) . it seems reasonable to 
apply our knowledge of Millennium projection effects to 
our SDSS pair statistics. Multiplying the SDSS redshift 
space N c (M r ) by the Millennium fzD{M r ) yields our best 
estimate of SDSS close pair statistics in real space. These 
results are given in Figure [6l as well as in Table [3] We 
find that the SDSS real space N c (M r ) is approximately 
0.011 for -22 < M r < -20, and then drops to about 
0.0055 for -20 < M r < -18. Given the larger error bars 
in our SDSS measurements, we cannot state with any 
certainty where the real-space N c (M r ) peaks; however, 
our results are consistent with a peak at M r < —20, as 
seen in the Millennium simulations (§ 15. ip . 

5.3. Comparison with Other Studies 

Most observed samples of close pairs have been too 
small to permit anything meaningful to be learned about 
the luminosity (or mass) dependen ce of close pair statis- 
tics. The most notable exception is lXu et all (|2004l ). who 
measure the pair fraction as a function of i^" s -band abso- 
lute magnitude for galaxies in the 2MASS survey. While 
they generally probe higher luminosities (—24.5 < Mk < 
—22.5) than we do, their two faintest bins (centered on 
Mk — —22) correspond roughly to our brightest SDSS 
bin (M r = —21.5); in this regime, they find a pair frac- 
tion of about 0.011 ± 0.005, which is consistent (within 
la) of our N c ~ 0.019±0.002 (Table©, despite some dif- 
ferences between their pair criteria and ours. However, 
IXu et all (|2004h find that the pair fraction increases to- 
wards higher luminosities (albeit with large error bars), 
whereas we see a hinted of a decrease with luminosity in 
this regime with Millennium. On the other hand, our ob- 
se rved trend is consistent wit h the semi-analytic results 
of iKhochfar & Burkertl (|2001h . who find that the z = 
merger fraction decreases as mass increases, over a range 
in masses which corresponds to M r < M*. 

We also compare our results with those of Berrier et al.l 
(2006) , who measure N c using a combination of N-body 
simulations and semi-analytic models. They find that 
N c increases steadily as the cumulative co-moving num- 
ber density increases. This is equivalent to an increase in 
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N c towards fainter limiting M r , over a range correspond- 
ing to M r > M* (us ing the SDSS LF measurements of 
iBlanton et alj (|2003aD l. We instead see a steady decrease 
in N c towards fain ter M r over this range. However, 
iBerrier et all (|2006f > include all companions brighter than 
the given number density (or equivalently, M r ). With 
this approach, the cumulative number of close compan- 
ions for a given galaxy can only go up as M r becomes 
fainter. It is therefore not surprising that their mea- 
surements of N c rise towards fainter M r . Given that we 
require companions to be within a factor of two in lumi- 
nosity of their host galaxy, there is no reason to expect 
ag reement between ou r measurements of N c and those 
of IBerrier et al.l (|2006h . We note, however, that our def- 
inition of N c is likely to be a better tracer of the major 
merger rate. 

6. FROM CLOSE PAIR STATISTICS TO THE MERGER 
RATE 

6.1. The Number Density of Close Companions 

We have measured the number of companions per 
galaxy (N c ), which tells us which galaxies are most likely 
to have close companions. The real-space distributions 
seen in Figure [6] indicate that low luminosity galaxies are 
the least likely to have close companions, and presumably 
the least likely to undergo imminent mergers. However, 
low luminosity galaxies are much more common than lu- 
minous galaxies, and therefore it is still possible that 
most major mergers may occur between low luminosity 
galaxies. 

We investigate this question by measuring the number 
of close companions per unit co-moving volume, here- 
after n c . In order to compute this quantity, we multiply 
the number of close companions per galaxy (N c (M r ) by 
the co-moving number density of galaxies (n(M r )). For 
SPSS, we estimate n( M r ) using the LF measurements of 
IBlanton et al.l (|2003af ). converting from the 01 r filter to 
the r filter using their recommended M r = Mo.i r — 0.16. 
For Millennium, we measure the galaxy number den- 
sity dire ctly from the ICroton et alj (l2006f) galaxy cata- 
log^UkeEiirin (|2007ir we find that the IBlanton et al.l 
(|2003a| ) LF matches the ICroton et ail d2006| catalog 
quite well overall, though the simulations overpredict the 
number of galaxies at fainter luminosities (the excess is 
~ 50% at M r — —18). We note that this excess helps to 
explain why N c (M r ) for wide separation pairs is steeper 
for Millennium than for SDSS (see SIQj). 

The resulting measurements of n c {M r ) are given in 
Figure [7] and in Tables [3] and [U Excellent agreement 
is seen between SDSS and Millennium. When summed 
over the range -22 < M r < -18, the SDSS n c = 
(2.10 ±0.13) x 10 _4 /i 3 Mpc~ 3 , which is consistent (within 
one sigma) with the Millennium n c = (1.98 ± 0.02) x 
10~ 4 ft. 3 Mpc~ 3 . In addition, the shapes of the distribu- 
tions are very similar, with n c ~6x 10 -5 /i 3 Mpc~ 3 mag -1 
for M r > —21, and a sharp drop in n c at higher lumi- 
nosities. For M r > —21, it therefore appears that the 
decrease in the real-space N c (M r ) towards fainter abso- 
lute magnitudes is balanced by a corresponding increase 
in host galaxy number density. On the other hand, the 
relative scarcity of luminous host galaxies leads to the 
rapid decline in n c towards bright absolute magnitudes. 
From Figure [7l we conclude that at least 90% of all ma- 



jor mergers occur between galaxies which are fainter than 
M r = —21. We note that, while there are no other di- 
rectly comparable measurements in the literature, the 
shape of our observed n c {M r ) is qualitatively similar 
to the shape of th e me rging galaxy mass fu nctions of 
iBundv et al] (|2005h and lHopkins et all (|2006h . 

6.2. The Galaxy Merger Rate 

We have reported estimates of the number of compan- 
ions per galaxy (iV c ) and per unit volume (n c ). With ad- 
ditional assumptions, both can be converted to merge r 
rates. Following the formalism of iPatton et alj {2000), 
we will refer to these as the galaxy merger rate (here- 
after TZ mg ) and the volume merger rate (7£ mg ) respec- 
tively. Assuming that two real-space close companions 
(ie., one galaxy pair) lead to one merger, and that the 
average timescale for such mergers is T mg , it follows that 
Ti-mg = 0.5-/V c /T m g and 7l mgv = 0.5n c /T mg , where both 
N c and n c refer to real-space measurements (Figures [6] 
andEJ. 

To proceed further requires an estimate of T mg . We be- 
gin by taking the simplest approach, which is to assume 
that T mg is equal to the dynamical friction timescale 
of a typical pair in the sample (after projection effects 
have bee n removed). We ad opt the estimate of 0.5 Gyr 
given bv IPatton et~aT1 (|2000f l. since our pairs are chosen 
with similar criteria. This timescale estimate is com- 
parable to or somewhat larger tha n several more re- 
cent estimates given in the literature (Ivan Dokkum 2005; 
iBell et al.1 l2006at iDe Propris et alJ l2007h . While the 
merger timescale may depend on mass (and hence abso- 
lute magnitude), N-body simulations appear to indicate 
that mass ratio is the most significant factor (jJiang et al.l 
2008) , and this should be independent of M r for the pair 
criteria we employ. 

The resulting merger rates are plotted using the right 
hand axes in Figures [6] and [7] The galaxy merger rate 
is found to peak at about 0.01 mergers per galaxy per 
Gyr, while the volume merger rate plateaus at about 
6 x 10~ 5 /i 3 Mpc Gyr mag -1 . As the conversion from 
pair statistics (iV c and n c ) to merger rates (72. mg and 
K mgv ) is independent of M r , the trends with M r de- 
scribed in Sections [5] and 16.11 apply here as well. We 
again conclude that while galaxies brighter than L* have 
the highest likelihood of being involved in major mergers, 
most major mergers take place between galaxies fainter 
than L*. 

How do these merger r ate estimates compare with 
others in the literature? IDe Propris et al.1 (|2007h es- 
timated the merger rate using two methods: close 
galaxy pairs (where both galaxies were required to have 
—21 < Mb < —18) and high asymmetries (as an in- 
dicator of ongoing mergers). They found a merger 
fraction close to 2% for both methods, which corre- 
sponds to 72. m g ~ 0.02 Gyr -1 using our T mg = 0.5 
Gyr. This is approxim ately twice as high as our result. 
IDe Propris et a l. (2007) also report a volume merger rate 
of (5.2 ± 1.0) x 10- 4 /i 3 Mpc -3 Gyr _1 , whereas we find 
TZ mgv = (1.4 ±0.1) x 10- 4 /i 3 Mpc^ 3 Gyr _1 when sum- 
ming over a comparable range in absolute magnitude 
(—22 < M r < —19). After recomputing their result with 
our T mg = 0.5 Gyr (instead of their 0.3 Gyr), we again 
find that their merger rate is roughly double ours. The 
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most likely explanation for this d ifference is that both the 
pairs and asymmetry methods of lDe Propris et all (|2007l ) 
include both major and minor mergers, whereas our mea- 
s urements app l y to m ajor merger candidates only. 

iMaller et alJ ((20061 ) measure the galaxy merger rate 
as a function of mass, mass ratio, and redshift, using 
a smoothed particle hydrodynamics (SPH) cosmological 
simulation. For mergers above a 1:2 mass ratio (compa- 
rable to our 1:2 luminosity ratio), they find the galaxy 
merger rate at z — 0.3 to be lZ mg = 0.054 Gyr -1 for 
high mass galaxies. This range in mass corresponds ap- 
proximately to —18.8 < M r < —17.7 (derived from Ta- 
ble 1 of Mailer et al. (2006)). While they do not re- 
port precise values of lZ mg at lower redshifts, they do 
find that 7?. mg declines quickly towards z = 0, yielding 
T^mg i$ 0.02 Gyr -1 at z < 0.1. This value is consis- 
tent with the 7Z mg ~ 0.005 Gyr -1 we find at z ~ 0.03 
over this range in M r , part i cularl y given that no er- 
ror bars on the IMaller et al.l (|2006[ ) estimate are avail- 
able. Interestingly, they find that the merger rate is 
lower for lower mass galaxies, in agreement with the 
trends we find in F igure [6] (particularly for Millennium). 
IMaller et all {2006) also measure a volume merger rate of 
about 1 x 10~ 4 /i 3 Mpc~ 3 Gyr _1 (with large error bars) for 
high mass galaxies at z — 0.1; this is consistent with the 
TZ mgv ~ 6 x lQ-^Mpc^G yr" 1 we find in FigureEl In 
general, IMaller et al.l ((2006) find that 1Z mgv declines to- 
wards lower masses (though this effect is significant only 
at z > 0.4), implying that one might expect to see 7£ rog 
decline at fainter absolute magnitudes than probed by 
o ur SPSS and Millenn ium catalogs. 

iMasiedi et alJ (|2006l ) measured the volume merger rate 
for pairs of luminous red galaxies (LRG's) in the SDSS, 
and found K mgv < 0.6 x 10 4 Gyr" 1 Gpc -3 . Both mem- 
bers of their LRG's pairs were required to have abso- 
lute magnitudes of —23.2 < M g < —21.2, which is dif- 
ferent in nature from our major merger requirement of 
|AM r | < 0.753. While our SPS S pairs are less luminous 
than the IMasiedi et"aT1 (|2006l ) LRG's, our Millennium 
volume merger rate estimates in Figure [7] yield a similar 
result {K mgv < 1 x 10 4 Gyr" 1 /i 3 Mpc~ 3 ) i n the LRG lu- 
minos ity range. Finally, we note that while Masiedi et al. 
(2008) probe the luminosity dependence of close compan- 
ions of LRG's in the SPSS, their (LRG) host galaxies are 
more luminous than the host galaxies we probe here, and 
most of their close companions are not luminous enough 
to qualify as major merger candidates; therefore, there 
is essentially no overlap between our studies. However, 
their Figure 3 is consistent with a major merger rate of 
TZ mg ~ O.OOlGyi'- 1 for typical LRG's, which is in quali- 
tative agreement with the hinted drop of N c (M r ) towards 
high luminosities which we describe in § 15.31 

We caution the reader that our estimates of the merger 
rate depend on two quantities: the close pair statistics 
we have measured (N c or n c ), and the merger timescale 
(T mg ) that we have estimated. The latter is considerably 
less certain, and different choices of the merger timescale 
and its dependence on mass (or luminosity) may signifi- 
cantly change our estimated merger rates. The most use- 
ful che ck on our results comes from lKitzbichler fc~ Whitc 
(2008), who use the Millennium simulation to devise fit- 
ting formulas for estimating the merger timescale for 
pairs within a given projected physical separation (here- 



after T merg0 ). While none of their fitting formulas are di- 
rectly applicable to our sample, the closest match comes 
from using their equation 9 with r p < 30 h^ 1 kpc, 
v p < 300 km s _1 , and z — 0. We approximate th e limit- 
ing stellar mass as the median lCroton et alJ (|2006f ) stellar 
mass of galaxies corresponding to the faintest allowable 
co mpanions for a given host g alaxy. We note also that 
the lKitzbichler fe White] (|2008l ) fitting formulas apply to 
pairs in which galaxy stellar masses differ by at most a 
factor of 4, whereas our sample consists of pairs in which 
galaxy luminosities differ by a most a factor of 2. 

These choices yield T mevge estimates which increase 
monotonically from about 1 Gyr at M r ~ —23 to 4.1 
Gyr at M r ~ —17.5 (the range covered by our Millen- 
nium sample). At face value, these results may appear 
to be at odds with the fixed merging timescale of 0.5 Gyr 
that we adopted earlier in this sec tion. However, there 
is an important difference: the iKitzbichler fc White 
(2008) timescale is relevant for projected pairs, while our 
timescale is for real space pairs. To make a direct com- 
parison, our timescale of 0.5 Gyr needs to be divided by 
fm(M r ); this leads to values of T merge ranging from 0.8 
Gyr to 3.2 Gyr over the same range in M r , which compare 
quite favorably with the preceding calculations. There- 
fore, de spite significant mismatches between our sample 
and the lKitzbichler fc White! (|2008f ) fitting formula crite- 
ria, both approaches lead to a similar luminosity depen- 
dence of the merging timescales for projected pairs, and 
the overall timescales are comparable. Moreover, our 
estimates of fso(M r ) (see Figure [5]) provide additional 
insight into the nature of the increase in T merge with de- 
creasing luminosity; namely, that this trend may be a 
consequence of most lower luminosity pairs having sepa- 
rations which are too large for merging to take place, 
rather than such systems simply undergoing a slower 
merging process. Nevertheless, it is clear that more work 
is needed in order to more accurately model the merging 
timescales of observed pair samples. 

6.3. The Remnant Fraction 

Given the above measurements of the merger rate, it 
is in principle possible to asses s the merging histor y of 
low redshift galaxies. Following iPatton et al.l (|2000f ). we 
estimate the fraction of galaxies which have undergone 
major mergers since z = 1. This quantity, which is re- 
ferred to as the remnant fraction (/ rom ), is given by 



/ rem (M r ) = 1 



N 1 



N c (M r ,Zj) 



A = l l-0.5A c (M r ,^)' 



(4) 



where N c (M r ) is measured in real space (e.g., §0, Zj 
corresponds to a lookback time of jT mg , and N refers 
to the number of merger timescales over the range < 
z < 1 (N = 15 for our chosen cosmology). In principle, 
one needs to measure N c (M r ) at < z < 1 in order to 
accurately compute f rcm (M r ). As we have only our low 
redshift (z < 0.1) measurements of N c (M r ) to work with, 
we make the simplifying assumption that N c (M r ) does 
not evolve with redshift. If instead N c (M r ) rises with 
redshift, as indicated by numerous studies (see jjU), the 
resul ting remnant fractions will be higher ((Patton et al.l 
2002). And of course, if the shape of N c (M r ) evolves 
strongly, the shape of our estimated f Tem {M r ) will also 
be in error. 



The Luminosity Dependence of the Galaxy Merger Rate 



9 



Nevertheless, we present our no-evolution estimate of 
the luminosity dependent remnant fraction in Figure [8] 
and in Tables [3] and [U The SDSS remnant fraction peaks 
at 9% for M r — -20.5 (roughly M*), and decreases to 
4% over the range —20 < M r < —18. Similar trends are 
seen for Millennium, although the Millennium remnant 
fraction drops to ~ 2% at M r = —18. 

6.4. The Luminosity Density of Close Companions 

We demonstrated in Section 16.11 that our observed 
n c (M r ) implies that most major mergers occur between 
galaxies which are fainter than L*. Of course, a ma- 
jor merger between two intrinsically faint galaxies affects 
much less stellar mass than a major merger of two lumi- 
nous galaxies. Therefore, in terms of the overall luminos- 
ity density of stellar mass that is participating in mergers, 
the number density of close companions does not tell the 
full story. To address this question, we instead measure 
the luminosity of close companions per unit co-moving 
volume (hereafter l c (M r )), using luminosity as a proxy 
for stellar mass. We compute l c {M r ) as follows: 

l c (M r ) = L c (M r )f 3d (M r )n(M r ), (5) 
where L r ( M r ) is the luminos ity in close companions 
per galaxy (jPatton et al.ll2000l ) as measured in redshift 
space, /3d(M r ) is as described in § 14.31 and n(M r ) is 
the co- moving number density of galaxies (see § 16. ip . 
When calculating l uminosities, we use M r (©) = 4.64 
(|Yasuda et al.ll2001[ ). 

The resulting measurements for both SDSS and Mil- 
lennium are given in Figure [9] and in Tables [3] and 
01 Excellent agreement is again seen between SDSS 
and Millennium. Both surveys exhibit a clear peak in 
l c (M r ) at (or very close to) M*, with a peak value of 
l c ~ 8 x 10 5 L Q /i 5 Mpc _3 mag _1 . The sharp decline to- 
wards the bright end is due to the lower number density 
of luminous host galaxies, while the decline towards the 
faint end is driven by the lower numbers and luminosities 
of companion galaxies. This distribution clearly implies 
that galaxies which have luminosities close to L* are the 
most relevant in terms of the overall involvement of stel- 
lar mass in major mergers. 

7. CONCLUSIONS 

We have measured the number of close companions per 
galaxy (N c ) as a function of absolute magnitude for both 
the SDSS and the Millennium simulation. For SDSS, 
we construct samples of host and companion galaxies, 
and correct for spectroscopic incompleteness. For Mil- 
lennium, we create a suite of mock redshift catalogs, av - 
eraging over different views of the ICroton et all §006) 
cube. Using close pair criteria designed to identify im- 
minent major mergers (5 < r p < 20 h^ 1 kpc, Av < 500 
km s _1 , and |AM r | < 0.753), we find general agreement 
between the observations and simulations. In redshift 
space, N c ~ 0.02 over the range -22 < M r < -18. The 
flatness of this distribution indicates that N c (M r ) does 
not simply trace the number density of galaxies; instead, 
small scale luminosity-dependent clustering appears to 
counteract this effect. 

Using three dimensional positions and velocities avail- 
able from the Millennium simulations, we measure the 
degree to which the detected galaxy pairs are contami- 
nated by projection effects, and find that the contami- 
nation is a strong function of absolute magnitude, rising 



from 45% at M r < M* to - 83% at M r = -18. We 
remove this contamination from both Millennium and 
SDSS pair statistics, yielding N c (M r ) measurements in 
real space. These measurements indicate that galaxies 
with M r < M* are the most likely to be undergoing 
major mergers at low redshift. However, by computing 
the number density of close companions (n c (M r )) in real 
space, we conclude that at least 90% of all major mergers 
occur between galaxies which are fainter than M* . With 
additional assumptions, we also estimate the galaxy and 
volume merger rates, which trace the real-space N c (M r ) 
and n c (M r ) respectively. We estimate that at least 8% of 
L* galaxies are likely to have undergone a major merger 
since z = 1, while this remnant fraction appears to be 
~ 4 times smaller for 0.1 L* galaxies. Finally, our mea- 
surements of the luminosity density of close companions 
indicate a clear peak at M*, implying that L* galaxies 
dominate in terms of the overall amount of stellar mass 
involved in major mergers at low redshift. 

Together, these results indicate that the low redshift 
merger rate depends strongly on luminosity (and pre- 
sumably mass) . This has a number of important impli- 
cations. For example, one would not expect the merger 
rates of massive galaxies (e.g., Masjedi et al. (2006)) to 
agree with those of L* galaxies (e.g., De Propris et al. 
2007). Also, the increase in projection effects for fainter 
galaxies indicates that samples of luminous galaxy pairs 
are more likely to provide bona fide merger candidates 
than samples of lower luminosity pairs. This is relevant 
if one wishes to assess the impact of merging on the con- 
stituent galaxies (e.g., enhanced star formation or dis- 
turbed morphologies). We also conclude that at low red- 
shift, recent merging history is likely to be most impor- 
tant for galaxies which are relatively luminous. Finally, 
given the clear peak in l c (M r ), it appears that galaxies 
which are much more or much less luminous than L* are 
unlikely to play an important role in the overall evolution 
of galaxies via major mergers. 
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TABLE 1 
SDSS Host Galaxies 



Host ObjlD 



T3S 



1T r 



Companion ObjlD 



587724198275842071 
587724197207670886 
587724197207670885 
587731513153159200 
587731513153159199 



10.610426 14.70509 0.05539 -19.72 587724198275842073 

23.645067 13.18706 0.06304 -20.17 587724197207670885 

23.646278 13.18500 0.06226 -20.50 587724197207670886 

42.731424 0.36690 0.04384 -20.10 587731513153159199 

42.739416 0.36940 0.04423 -19.76 587731513153159200 



TABLE 2 

SDSS Redshift Space N c (M r ) as a Function of Projected Pair Separation 

r v range -22 < M r < -21 -21 < M r < -20 -20 < M r < -19 -19 < M r < -18 

(h- 1 kpc) (z = 0.107) (z = 0.074) (z = 0.048) (z = 0.031) 

5-20 0.0187 ± 0.0018 0.0238 ± 0.0016 0.0172 ± 0.0018 0.0226 ± 0.0034 

20-40 0.0225 ± 0.0018 0.0275 ± 0.0016 0.0298 ± 0.0020 0.0284 ± 0.0026 

40-60 0.0174 ± 0.0015 0.0332 ± 0.0015 0.0345 ± 0.0016 0.0329 ± 0.0024 

60-80 0.0180 ± 0.0013 0.0346 ± 0.0011 0.0389 ± 0.0016 0.0398 ± 0.0026 

80-100 0.0210 ± 0.0011 0.0360 ± 0.0011 0.0405 ± 0.0016 0.0470 ± 0.0028 



TABLE 3 
SDSS Close Pair Statistics 



M r N c f 3D N c n c / rcm l c 

(redshift-space) (%) (real-space) (%) 

^2lM 0.0187 ± 0.0018 55.3 ± 0.4 0.0103 ± 0.0010 1.52 ± 0.15 7.5 ± 0.7 3.44 ± 0.37 

-20.50 0.0238 ± 0.0016 52.5 ± 0.3 0.0125 ± 0.0008 6.93 ± 0.44 9.0 ± 0.6 7.97 ± 0.58 

-19.50 0.0172 ± 0.0018 33.8 ± 0.6 0.0058 ± 0.0006 5.78 ± 0.60 4.3 ± 0.4 2.81 ± 0.33 

-18.50 0.0226 ± 0.0034 22.8 ± 0.7 0.0052 ± 0.0008 6.78 ± 1.04 3.8 ± 0.6 1.57 ± 0.26 



TABLE 4 
Millennium Close Pair Statistics 



Mr 


N c 

(redshift-space) 


(%) 








N c 

(real-space) 




n c 






/rem 

(%) 






lc 




-22.75 


0.0142 ± 0.0010 


45.9 ± 


3 


.0 





.0068 ± 0.0007 


0.03 


± 





.00 


5.0 ± 0. 


5 





.17 





01 


-22.25 


0.0173 ± 0.0004 


57.9 ± 


1 


.2 





.0101 ± 0.0003 


0.18 


± 





.01 


7.3 ± 0. 


.2 





.79 





02 


-21.75 


0.0163 ± 0.0002 


54.9 ± 





.5 





.0089 ± 0.0001 


0.71 


± 





.01 


6.5 ± 0. 


.1 


2 


.05 


0. 


.02 


-21.25 


0.0185 ± 0.0002 


55.5 ± 





.6 





.0103 ± 0.0001 


2.52 


± 





.03 


7.5 ± 0. 


1 


5 


.02 





.06 


-20.75 


0.0207 ± 0.0002 


55.5 ± 





.4 





.0115 ± 0.0001 


5.57 


± 





.06 


8.3 ± 0. 


1 
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.76 


0. 


.07 


-20.25 


0.0191 ± 0.0002 


48.5 ± 





.6 





.0093 ± 0.0002 


6.85 


± 





.12 
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.58 


0. 


07 
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.8 
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0. 
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.95 
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.9 
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.35 


2.2 ± 0. 


1 
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2 





.56 





.04 



12 



D. R. Patton & J. E. Atfield 




0.05 0.1 0.15 

z 

Fig. 1. — Absolute magnitude is plotted versus redshift for 10,000 galaxies selected at random from our SDSS spectroscopic catalog of 
about 337,000 galaxies. The dashed (black) lines enclose all potential companion galaxies; galaxies within this region have 14.5 < m r < 17.5. 
The solid (red) lines enclose all potential host galaxies; galaxies within this region differ by at least 0.753 magnitudes from the dashed lines. 
For these potential host galaxies, all companions within a luminosity ratio of 1:2 are therefore detectable. 
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Fig. 2. — The ratio of spectroscopic to photometric pairs (N zz /N pp ) is plotted versus angular separation (8), with error bars computed 
using the Jackknife technique, (a) The upper panel shows the data before corrective weights have been applied, and a substantial deficit 
of spectroscopic pairs is seen at small separations (8 < 55"). At larger separations, N zz /N pp converges to the square of the overall 
spectroscopic completeness of the sample (/^ ~ 0.78), as indicated by the dashed line. The solid lines show our model fit to the data, 
which is used to correct for the small scale spectroscopic incompleteness, (b) The lower panel shows the data after these weights have been 
applied. The corrected data points provide an excellent fit to the overall spectroscopic completeness (solid line), confirming that the small 
scale incompleteness has been successfully removed. 
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Fig. 3. — N c is plotted versus absolute magnitude for SDSS, for pairs in five ranges of projected separation. Over the range —22 < M r < 
— 18, the mean N c for close pairs (5-20 h~ kpc) is 0.021 ± 0.001. The N c distribution is approximately flat for the smallest separation 
pairs, and becomes steeper as pair separation increases. 
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Fig. 4. — N c is plotted versus absolute magnitude for the Millennium simulation, for pairs in five ranges of projected separation. As with 
SDSS, the N c distribution is found to be approximately fiat for the smallest pair separations. Over the range —22 < M r < —18, the mean 
N c for close pairs is 0.0183 ±0.0001, which is about 15% lower than found for SDSS. The Millennium N c distribution steepens towards the 
faint end as pair separation increases. This is similar to what was seen for SDSS, but is considerably steeper. 
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Fig. 5. — Using the Millennium simulation, the fraction of companions which are close in three dimensions (/3d) is plotted versus absolute 
magnitude for close pairs (5 — 20 h^ 1 kpc). Projection effects are seen to increase towards fainter absolute magnitudes. 
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Fig. 6. — N c (M r ) is plotted in both redshift space and real space, for close pairs (5 — 20 ft -1 kpc). For both SDSS and Millennium, 
conversion from redshift space to real space pair statistics (ie., the removal of projection effects) was carried out by multiplying N c (M r ) 
by the Millennium /3£i(M r ). The right hand axis displays the corresponding scale for the galaxy merger rate (7?. mg ), and is relevant only 
for measurements in real space. 
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Fig. 7. — The number of close companions per unit co-moving volume (n c ) is plotted versus M r for SDSS and for Millennium. In both 
cases, n c is computed by multiplying the real-space N c by the number density of galaxies. Excellent agreement is seen between the two 
samples. We conclude that most mergers occur between galaxies with M r > —21. The right hand axis displays the corresponding volume 
merger rate. 
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Fig. 8. — The remnant fraction /rem, expressed as a percentage, is plotted versus M r for SDSS and for Millennium. In both cases, 
the error bars are computed using only the uncertainties in N c . Given the assumptions that go into the calculation of /rem, the true 
uncertainties are likely to be considerably larger. 
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Fig. 9. — The luminosity in close companions per unit co-moving volume (Z c ) is plotted versus M r for SDSS and for Millennium. Excellent 
agreement is seen between the two samples. l c (M r ) exhibits a clear peak at or very close to M* . The decline towards brighter absolute 
magnitudes results primarily from the low number density of such systems, whereas the decline at the faint end results from the decreasing 
luminosity of companions. In terms of stellar mass (as inferred from galaxy luminosities), this implies that L* is the most relevant regime 
for major mergers. 



